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Abstract 

i-^ \ A classic "no-go" theorem in one-dimensional quantum mechanics can be evaded when the 

potentials are unbounded below, thus allowing for novel parity-paired degenerate energy bound 
states. We numerically determine the spectrum of one such potential and study the parametric 
variation of the transition wavelength between a bound state lying inside the valley of the potential 
and another, von Neumann- Wigner-like state, appearing above the potential maximum. We then 
construct a modified potential which is bounded below except when a parameter is tuned to vanish. 
We show how the spacing between certain energy levels gradually decrease as we tune the parameter 
to approach the value for which unboundedness arises, thus quantitatively linking the closeness of 
degeneracy to the steepness of the potential. Our results are generic to a large class of such 
potentials. Apart from their conceptual interest, such potentials might be realisable in mesoscopic 
systems thus allowing for the experimental study of the novel states. The numerical spectrum in 
this study is determined using the asymptotic iteration method which we briefly review. 
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I. INTRODUCTION 



Bound states in quantum mechanical potential problems usually have energies below the 
potential maximum, as most quantum mechanics textbooks remind us. However, already in 
1929 von Neumann and Wigner l| first noted the theoretical existence of a class potentials 
which support normalizable bound states above the potential maximum, embedded in the 
continuum of scattering states. Later work by Stillinger and Herrick [2j provided a more 
detailed understanding of such unusual bound states, indicating their essential quantum 
mechanical nature. The von Neumann-Wigner above-barrier bound states, were eventually 
experimentally verified in the 1990s, by Capasso et al. [3], through observations on allowed 
transitions, between a state inside the valley and another above the potential maximum, in 
semiconductor heterostructures. 

Recently, potentials unbounded from below, in one dimension, have attracted interest 
because of certain unusual features in their energy eigenvalues and eigenfunctions. The class 
of potentials studied were non-singular in any finite domain but asymptotically approached 
negative infinity. It was shown that such potentials not only have von Neumann-Wigner 
states, they also exhibit a new feature: parity-paired degenerate bound states with real 
energies {4, 5], a property prohibited for regular one-dimensional potentials [6]. The study 
of such potentials is thus conceptually interesting as they highlight how commonly accepted 
"no-go" theorems may be evaded. 

It should also be noted that such unbounded-below potentials are not just mathematical 
constructs, but have appeared in the context of localization of fields on the 3-brane in 
the study of the so-called braneworld models with warped extra dimensions in high energy 
physicsj?]. As mentioned in J4], and discussed below, such potentials might be realisable, to 
any given degree of accuracy, in actual mesoscopic systems. 

In this paper, we study in detail the spectrum of two potentials so as to anticipate 
future experimental investigations as suggested in Ref . j4j . We first numerically compute the 
eigenvalues, and hence transition wavelengths, for the cosh-sechpotential to supplement the 
partial information that is available in earlier analytical work (4]. 

We then modify the unbounded potential in such a way that the unboundedness appears 
only when we tune a parameter in the potential to vanish. This allows us to study how the 
degenerate states of the unbounded potential split as the parameter is varied, and as one 
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moves to the more realistic bounded potential. Since in [4j we had shown the existence of 
large classes of potentials supporting degenerate bound states, the splitting of the degenerate 
states for the bound versions of those potentials is naturally expected and thus our results 
are qualitatively generic. 

The eigenvalues are determined with the help of the asymptotic iteration method (AIM), 
a relatively novel technique to solve second-order linear homogeneous differential equations 
with variable coefficients js ij] (details of the method are provided in Appendix A). Its 
asset is that it is able to easily circumvent some of the profound analytical and numerical 
difficulties posed for the Schrodinger equation with potentials which are either singular or 
unbounded from below. 

Our article is organised as follows. In Section II, we discuss the cosh-sech potential and 
determine its spectrum using the asymptotic iteration method (AIM). In Section III we 
discuss the degeneracy of states and how the levels split in the modified bounded potential 
that we construct. Section IV summarises our results while the appendices (Appendix A 
and Appendix B) contain some technical details about AIM and physical units, respectively. 

II. A MODIFIED POSCH TELLER POTENTIAL 

Consider the cosh-sech potential, 



b 2 ( r 

V(x) = — — cosh 2 x — I a 2 — — ) sech 2 x (1) 



4 V 4 

This differs from the potential studied in [3| just by an additive constant, —b 2 /4, since they 
131 ] choose the first term as a sinh function. 
The above potential has a maximum at 



X r 



cosh 1 



4a 2 -IN 1/4 



b 2 



which exists only for b < ^Aa 2 — 1 since cosh(x) > 1. The extrema of the potential function 
are V max = — |V4a 2 — 1 and Vmin = — (o 2 — |) — ^. The potential is of the same form as 
the volcano potential of Ref . [5| , 



V(x) = -A x cosh 2 " x - A 2 sech 2 x (2) 



with v = 1, which can be generated from the general ansatz of Ref. {4] by choosing f(x) 
cosh""/ 2 x. 
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To solve the time-independent Schrodinger equation using the AIM (see Appendix A) 



the wave function is written as ip = e f^y 



131 ] . where 



f(x) = ib smh(x) + (a — 1/2) lncosh(x) 



(3) 



A further change of independent variable to u — sinhx, transforms the Schrodinger equation 
to 

y" = \ (u)y' + s (u)y (4) 

where 



A («) = 2 



ib 



a — 1 

u 2 + 1 



-u 



S {U 



2iub(a-l) E a- a 2 -1/4 3b 2 



u 2 + 1 



u + 1 



(5) 



(6) 



The above form of the Schrodinger equation is necessary in order to use AIM, to find energy 
eigenvalues for different a, b values. 



A. Energy spectra and transition wavelengths 

We now present the results of our calculations using the asymptotic iteration method. 
All computations were carried out using standard routines in Mathematica with default 
precision. 

In Fig. [H we plot the dimensionless energy eigenvalues (see Appendix B for a discussion 
on physical units) for the unbounded potential with a = 20 and b varying from up to 10. 
One can observe that for energy eigenvalues, as b increases, energies increase in general. Our 
calculations (not shown here) also show that the energy eigenvalues monotonically decrease 
with the value of the parameter a. 

The transition wavelength between a pair of states can be determined from the difference 
in their energies, after converting to suitable physical units, as outlined in Appendix B. 
Fig. [2] demonstrates a pair of eigenstates located above (El) and below (E2) the maximum 
in the potential well. The parameters a and b of the cosh-sech potential have been chosen 
so that the transition wavelength lies in the visible range of the electromagnetic spectrum. 

The transition wavelength can be varied over the range of interest by modifying the 
parameters a and b; from Fig. [21 it is evident that for a = 10 to 16, A is in the IR region, 
and greater values of a lead to transitions in the visible region. Similarly, Fig. H] shows a 
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FIG. 1: Variation of energy eigenvalues with b for the cosh-sech potential for a = 20. 




FIG. 2: Schematic of the cosh-sech potential and the transition wavelength between a pair 
of states El and E2 (in red) located outside and within the potential well (in blue) 
respectively, for a = 19, b = 2. The wavelength A = he/ (El — E2) = 0.569 fim lies in the 
visible range of the spectrum. The energy scale, the potential and the x-coordinate are 

plotted in atomic units. 
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FIG. 3: The variation of the transition wavelength between a pair of states outside and 
inside the cosh-sech potential well with the parameter a, for 6 = 2. The black dashed line 
represents the rough boundary between the visible and the infra-red regions of the optical 

spectrum. 



variation of the transition wavelength from the IR to the ultra-violet region by changing the 
parameter b. Note that since the energy differences increase with an increase in b, hence A 
moves from the longer IR wavelengths to the relatively short visible region. 

For illustration, the length scale parameter xq (see Appendix B) had been chosen to be 

o 

10 A (= 0.001 jum) in the above calculations. Xq can be modified suitably to obtain orders 
of magnitude variation in the energy eigenvalues (in physical units) and the corresponding 
transition wavelengths. 



III. DEGENERACY 



As mentioned earlier, a curious feature of the class of unbounded potentials being consid- 
ered here is the existence of degeneracy of the bound eigenstates despite the fact that we are 
dealing with one-dimensional quantum mechanical problems % |5( . Not all potentials un- 
bounded from below as x — > ±oo support degenerate bound states; for instance the quartic 
anharmonic oscillator with V(x) = x 2 — Ax 4 , A > is not known to exhibit degeneracy. 
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b 

FIG. 4: The variation of the transition wavelength between a pair of states outside and 
inside the potential well with the parameter b, for a = 16. The black dashed lines represent 
the rough boundary between the infra-red, the visible and the ultraviolet regions of the 

optical spectrum. 



Here we show, numerically, that the degeneracy exists in our potential for many excited 
states which lie above the previously analytically obtained pair of states (possibly ground 
states) in Ref.j^]. 

Consider the cosh-sech potential with a = 3,6 = 1. After 
12 iterations, the following eigenvalues are obtained using the AIM: 
-6.47301,-6.3402,-2.6222,-2.6058,0.9607,0.9624. It can be inferred that the en- 
ergy eigenvalues occur in closely separated pairs. Moreover, the separation between the 
two eigenvalues in a pair keeps on decreasing with an increase in the number of iterations, 
suggesting that they are ideally degenerate. The energy eigenvalues below 
V(x) at x = do not correspond to bound states. This is discussed in [L 
where the authors show how the below-minimum states are not linked to total-transmission 
(TT) modes, unlike genuine bound states which are always constructed out of the TT 
modes. Further, we have also seen that these below-minimum states do not occur in pairs, 
a feature different from the above-minimum bound states. Thus, in all our evaluations 
using AIM we have ignored such states. 



;he minimum m 
I (see also Il3j), 
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FIG. 5: The gradual modification of the cosh-sech potential V c (x) (Eq. [7]) with the 
parameter c for a = b = 1. The parameter c can be used to tune the unboundedness of the 

potential. 

Next, we consider a modified potential which diminishes to zero instead of going to infinity 
±oo, and which should be realisable in a laboratory. We multiply the cosh term 
in EqJTJ which is responsible for the unboundedness, by an exponentially decaying term 
controlled by the parameter c. 

b 2 ( l\ 

V c {x) = — — cosh 2 x x exp(— cx 2 ) — I a 2 — - j sech 2 x (7) 

The original potential can be recovered by setting c = 0. As c increases, the first term gets 
suppressed, and the potential retreats to zero at a faster rate. This effect is illustrated in the 
Fig. E]for a = b = 1. For c = 0.20, the behaviour of V c (x) is quite close to that for c = near 
the origin, but far from x = 0, the potential decays to zero instead of the latter case which 
goes to — oo. This offers promise in investigating the relationship between unboundedness 
and degeneracy. On the other hand, the c = 0.5 case is similar to the bounded potential 
well or the Posch- Teller potential rather than the unbounded cosh-sech potential. Thus, c 
can be used to tune the steepness of the potential away from the origin. 

The energy spectrum of the modified potential V c (x) can be obtained via the AIM by 
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changing the expression for s (u) in Eq. 0to: 

. . 2iub(a — l) E a -a 2 -1/4 b 2 . . . , 2 , n . 

s o («) = T~T^ ~ "TXT + 211 + T 2 + exp _c arcsmh M 8 

U z + 1 M z + 1 IT + 1 4 

The arcsinh(w) factor occurs in the exponent due to the intermediate variable transformation 
u = sinhx. 

Fig. |6] shows the energy eigenvalues of the pair of degenerate states with the parameter 
c for the specific case of the potential considered above (viz. a = b = 1). As we have 
seen before, a bound state exists for the unbounded potential at -0.25. It is evident from 
the plot that the pair of degenerate levels splits into separate levels with an increase in c. 
The splitting is also shown separately in Fig. [7] where the energy difference is plotted. We 
note that as we move towards the left (i.e. c — > 0), the minute difference between the two 
degenerate eigenvalues which still remain even at c = 0, is due to numerical error in the AIM 
rather than actual splitting. This is why the curve is steadily increasing for c > 5 x 1CT 3 
and flattens for c < 5 x 1CT 3 . In fact, E\p — £a,i does not go below ~ 5 x 10~ 6 even for 
c = 0, if we use 12 iterations. 




0.1 0.2 0.3 0.4 0.5 

c 

FIG. 6: Splitting of degeneracy with an increase in the parameter c. As the potential 
varies, the difference in energy between the two eigenstates increases. 

Thus, tuning the parameter c modifies the potential from an unbounded to a bounded 
one, and gradually increases the separation between the levels which are degenerate in 
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FIG. 7: Variation of the difference AE as a function of c. The axes are logarithmic. 



the unbounded case. While it might be impossible to physically realize the unbounded 
potential (c = 0), its variants such as the potentials V c (x) with c > can be experimentally 
fabricated through methods of band-gap engineering in semiconductor heterostructures 15) , 
where potentials with multiple wells and barriers do arise and are not very uncommon. 
The experimental energy eigenvalues can thus be compared with the theoretically predicted 
results with the help of measured transition wavelengths. 



IV. CONCLUSION 



In this paper we first calculated the energy spectra and transition wavelengths between 
states with energies lying within the potential valley and above the potential maximum for 
the volcano potential (pQ), thus extending in detail the analytically obtained results in J4]. 

We emphasize that we have been able to demonstrate, through our numerical work in this 
paper, that there exist several degenerate excited states in the full spectrum, beyond the sin- 
gle pair obtained analytically in Ref.Q]. We then constructed a new potential, depending on 
a tunable parameter c, which was similar to the volcano potential near the origin but which 
was bounded as x — > ±00 for c ^ 0. This non-singular potential exhibited non-degenerate 
energy eigenvalues. For c = one has the original volcano potential with degenerate states 
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and as c is increased, creating a bounded potential, the splitting of degenerate states also 
increased. 

As stated in Section I, the properties of the bounded potential we studied here are likely 
to be qualitatively similar to others that can be constructed from the general ansatz of 
Ref. 4( and plausibly amenable to experimental realization in semiconductor heterostruc- 
tures. We believe that laboratory fabrication of systems where such potentials may arise 
could be a route to the observation of the novel parity paired (almost) degenerate states in 
such one-dimensional problems. Further, the splitting of the degeneracy will lead to tun- 
able (varying c) long wavelength transitions between the closely spaced split levels for the 
bounded potential. As for the case of the von Neumann- Wigner states, we feel that the 
study of such parity-paired degenerate states is of intrinsic conceptual interest. 

Finally, the present study once again shows that the AIM is a useful semi-analytical tool 
to solve Scrodinger-type eigenvalue problems. We have verified the pre-existing bound state 
analytical solutions for unbounded potentials, with the help of the AIM, and extended the 
method to determine the entire bound state energy spectrum of these potentials which are 
unbounded below. 



Appendix A: The asymptotic iteration method 

The asymptotic iteration method (AIM) [8] is a general semi-analytical technique for 
solving second order, linear homogeneous ordinary differential equations (ODEs). It has 
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been used quite extensively over the last few years, in a variety of contexts (see 
a few references). The basic idea behind the method is as follows. Let Ao(x) and sq(x) be 
functions defined over a certain interval, having appropriately many successive derivatives. 
The differential equation we are concerned with is given as 

% = ^) d f x + s {x)y (9) 

This equation is to be solved by the AIM. The technique consists of iteratively differentiating 
the equation to obtain an ODE of the same general type, where the coefficients of the next 
step Afc(x), Sfc(x) are related to the previous coefficients by the expressions 

Afc(x) = ^A fe -i(a;) + s k -i(x) + A (x)A fc -i(x) (10) 
Sk(x) = — s fc _i(x) + s (x)A fc -i(x). (k G IN) (11) 
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The iteration is terminated when the ratio between the coefficients becomes independent of 
the index k, 

TTT = \ TT = a W 12 

Afc(a;) Afc_i(a;) 

in which case the equation is exactly solvable with the help of the AIM and exact eigenvalues 
are obtained by finding the roots of the discriminant equation 

= s fc _i(x)A fc (x) - s fc (x)A fc _i(x) = (13) 

On the other hand, if the ratio of succesive coefficients is not strictly independent of k, exact 
solutions cannot be found by the AIM. Nonetheless, an approximation to the eigenvalues 
can be obtained by forced imposition of the condition in Eq. [12] for a sufficiently large value 
of k, wherein lies the asymptotic nature of the method. 

n 

Furthermore, after some algebra, it is shown in |8j that the solution y(x) to Eq. 0can be 
obtained by the following integral relation: 

ft 



y[x) 



exp f- J adt^j C 2 + Ci J exp (J (A (r) + 2a(r)) drj dt 



(14) 



The time independent Schrodinger equation is of the form 

-^ + —(E-V(x)^ = (15) 

By a suitable transformation of ip(x) which reflects the asymptotic behaviour of the wave 
function, Eq. [TBI can be recast into a form resembling Eq. [HI where Xq(x) and Sq(x) depend on 
the energy E. Applying the asymptotic iteration method, the A^'s and s^'s are determined 
to arrive at the "quantization" condition of Eq. [13j Since 5k depends on both x and E, 
it becomes necessary to choose a suitable value of x = xq in order to solve for the energy 
eigenvalue E. For exactly solvable models, an arbitrary choice of xq suffices to arrive at the 
correct eigenvalue. On the contrary, for models which are not exactly solvable, the choice 
of Xq affects the rate of convergence of the method, and hence the accuracy of the energy 
calculated. 

It is noteworthy that the AIM is essentially an analytical technique, as the successive 
derivatives required are determined by analytic or symbolic differentiation. However, the 
form of 5k(x) is usually so complicated that the use of numerical algorithms is necessary 
in order to calculate the roots of the equation 5 fc = 0, with a choice of x = x where x is 
usually taken to be the critical (maximum/minimum) point of the potential. 
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We tested the AIM on the potentials mentioned in J4] and verified that the numerical 
eigenvalues agree with the partial information already available analytically. For example, 
we have checked that for a = 1 there is always one energy eigenvalue E = —0.25 V6, which 
corresponds to the known exact solution J] . New results are presented in the main body of 
this paper. 



Appendix B: Physical units 

In dimensionless (atomic) units, the time independent Schrodinger equation is: 

d 2 



dx 1 



+ V{x) \ i){x) = E n ijj(x) (16) 



where x is a dimensionless variable proportional to the one-dimensional spatial coordinate 
and V(x) is the potential proportional to that in the actual quantum mechanical problem. 
On the other hand, the Schrodinger equation in real units is: 

h 2 d 2 



2/;jf// . 2 +V(r)|V(r) = «(r) (17) 

The dimensionless coordinate x can be related to the dimensional coordinate r with the help 
of a length scale parameter xq, x = t/xq. Here, Xq is a measure of the width of the potential. 
Substituting this relation in Eq. [T71 we get: 

-^2+ ^r^(xMr) = _^°^(r) (18) 

Comparing with Eq. [161 

2 -^V(r) = V(x), 2 -^E n = E n (19) 

Thus, we have, 

h 2 

E^real units) = ^E n (no dimensions) (20) 



2mx, 



For example, an electron with a mass m e = 9.11 x 10 kg has an energy: 

~ . , , 6.1042 x 10~ 39 £ n (no dim.) 
i? n (real units) = ^ J 

x 

3.8104 eV ^ . v 
= — -B„(nodam.) (21) 

xl (in A) 
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If we take x = 10 A, the dimensionless energy is related to the real energy by the factor, 

£ n (real units) = 38.104 meV x E n (= 6.1042 x 10~ 21 J x E n ) (22) 

Corresponding to an energy difference AE, the transition wavelength A (in /jm) = he/ AE 
1.2424/AE(in eV). As an example, consider the red and violet ends of the visible spectrum 
corresponding to wavelengths of 0.8 [im and 0.4 /im, which translate to dimensionless energy 
differences of 81.51 and 40.75 respectively. By suitable choice of the parameters in the 
potentials, it is possible to obtain transition wavelengths in the optical range or beyond. 
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